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The behavior of perturbations is studied in cosmological models which consist of two different 
homogeneous regions connected in a spherical shell boundary. The junction conditions for the 
metric perturbations and the displacements of the shell boundary are analyzed and the surface 
densities of the perturbed energy and momentum in the shell are derived, using Mukohyama's 
gauge-invariant formalism and the Israel discontinuity condition. In both homogeneous regions, the 
perturbations of scalar, vector and tensor types are expanded using the 3-dimensional harmonic 
functions, but the mode coupling among them is caused in the shell by the inhomogeneity. By 
treating the perturbations with odd and even parities separately, it is found, however, that we can 
have consistent displacements and surface densities for given metric perturbations. 
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^ , PACS numbers: 98.80.Cq, 98.70.Vc, 95.35.+d 

o : 

^ , I. INTRODUCTION 

c . 

' In the real Universe, it is known at present that there are many superclusters and void-like objects with 
00 ', nonlinear overdense and underdense regions which are distributed everywhere on various scales. Such objects 
may have complicated structures in general, but their evolution has often been studied in comparatively large- 
scale cases, using a simplified assumption of spherical symmetry. 

As general-relativistic examples of spherical treatments we have the analysis due to the Tolman-Bondi- 
Lemaitre inhomogeneous dust model [1, 2], a single-shell model consisting of two different homogeneous regions 
, connected with a spherical shell [3, 4], a self-similar model being described using the self-similar solutions [5, 6], 
QQ a model consisting of two homogeneous regions connected with the self-similar intermediate region [7, 8, 9], and 
so on. 

In order to study the dynamical behavior of the large-scale superclusters and void-like objects and their 
influence on the CMB anisotropy, moreover, it is necessary to consider the gravitational instability of these 
models. The general gauge-invariant formalism for linear perturbations in spherical symmetric inhomogeneous 
Qh' models was derived by Gerlach and Sengupta [10, 11, 12, 13] in four-dimensional spacetimes. Recently the 
Q junction of perturbations in two homogeneous regions was treated by Mukohyama[14, 15] and Kodama et 
5^ al.[16] in higher dimensional spacetimes on the basis of the Israel discontinuity condition [17]. In this paper, 
\ we study the behavior of linear perturbations in the four-dimensional cosmological models with a single shell, 
d • using Mukohyama's formalism. Here only one of his doubly gauge-invariant variables are used, and in both 
^ ' homogeneous regions the perturbations are classified into scalar, vector and tensor types in three dimensional 
• : space and expanded in terms of three dimensional harmonic functions [18, 19] in order to clarify the physical 
^\ image of perturbations, though he classified them in a space of the shell surface. Among them the mode coupling 
is caused by the inhomogeneity in the shell. By treating the perturbations with odd and even parities separately, 
it is found however that the consistent expressions of the shell displacement and the energy-momentum tensor 
in the shell can be derived for given metric perturbations in both regions. 

In §11, we describe briefly the single-shell model as the background model, in which the velocity of the 
shell is determined through the junction condition. In §111, we show the basic equations for gauge-invariant 
perturbations and the junction condition based on Mukohyama's formalism [14, 15]. In §IV, we derive the 
intrinsic perturbations of metric and extrinsic curvature in the shell, and show the displacements of the shell 
and the intrinsic energy-momentum tensor in the shell which are derived from the junction conditions, by 
treating the perturbations with odd-parity and even-parity separately and imposing a localization condition on 
the wave-number dependence of the shell displacement. §V is given to the concluding remarks. In Appendix 
A, the harmonic functions in a homogeneous background model are shown, which are necessary for the analysis 
of the description of junction conditions. In Appendices B ad C, auxiliary quantities J\^u and Q^i, for the 
perturbed gravitational flcld and the intrinsic energy-momentum tensor in the shell, respectively, are shown. 
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II. BACKGROUND MODELS AND THE JUNCTION CONDITIONS 

The background universe is assumed to consist of two spatially homogeneous regions connected by a spherical 
shell. The four- dimensional line-elements in the inner region V~ and the outer region V+ are described as 

ds^ = g%{dx>')±{dx^)± = [a±{ri±)f[-{dn±f + dll], (1) 

dll = jtj{dx%{dx^)± = {dx±f + [<J±{x±)?dn\ (2) 

where the suffices /i, y and i,j run from to 3 and from 1 to 3, respectively, x\. = X±) f±(x±) = sinxj., x-t 
and sinhx± for the spatial curvature K± = 1,0,-1, respectively, and dO? = dO"^ + siv? Od^p^ . The shell is a 
time-like hypersurface S given as 

x'^ = Z^{y^), (3) 

where y*^ (M = 0, 2, 3) denote the intrinsic coordinates in the shell. In the following we specify the coordinates 
in the form 

xi=Zi{y^), xl = y\ (4) 

where the suffices A and a take the values (0, 1) and (2, 3), respectively. When we define the tangential triads 
^(M)± ~ ^' ^' ^(M)± ~ dx'^/dy^ , the intrinsic metric is given by 

Smat = ^(M)±^(N)±9%- (5) 
If we specify y° as the proper time r in the shell, the components of the triads can be expressed as 

e^i = (axV%o,aa;VayO,0,0) = F±(1,^;±,0,0), (6) 

e'^2)± ^ (0,0,1,0) and = (0,0,0,1), where F± = 7±/a±, v± = dxj./dxj. (= dxj./dri±) and j± = 

l/\/l — (w±)^. Here we take the units of the light velocity c = 1. The unit normal vector n^. of S is defined by 
the conditions n^±e^j^^^_|_ = (M = 0, 2, 3) and n^±n^ = 1, and the components can be expressed as 

< = (-^^±,l,0,0)7±/a±. (7) 
The line-element in the shell in terms of the intrinsic coordinates is 

ds"^ = quNdy^dy^ = -dr^ + R^dii^, (8) 

where y° = T,y'^ = e,y^ = <f, and R = a+cr(x+) = a_a{x-)- 

Since the intrinsic geometry is regular, the induced metric should be continuous and we have 

qMN+ = qMN- = QMN- (9) 

The extrinsic curvature of S is defined by 

Kmn± = ■^e.'^M)±^\N)±^ng^u = e('M)±e(W)± + '^'^■,1^]^ (10) 

where C represents the Lie derivative, n denotes the unit normal of S, and a suffix ; fj, denotes the four- 

dimcnsional covariant derivative with respect to .x^. For the energy- momentum tensor T'^" , the corresponding 
surface component S^,^ in the shell is defined using the Gaussian normal coordinates as 

S^^ = lim / T^^dC, (11) 

where C is the radial coordinate and C = is for the shell surface. The expression Smn in the intrinsic 
coordinates is given by 

Smn = (^^m)^\n)S iJ,y ^ (12) 
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where e^^-j is the triad in the gaussian normal coordinates. 

From the junction conditions derived by Israel, we obtain the relations between the jump of the extrinsic 
curvature and the surface energy- momentum tensor Smn as 

[Kmn]^ = -h^{Smn - qmn), (13) 



together with the additional relations 

'^M\N = L-'M 



and 

1 
2 

where = $+ - S = q'^^ Smn, is the inverse of the induced metric qmn, = SttG, and 



{K^^ + KfJ_)S^ = [T:]^, (15) 



T2i±=n^±elM)±{TI^)± T^^ = n^±nl{TI^)±, (16) 
Sm — Q^^Slm- (17) 

In the case of dust matter, we have 

T^iv± = P±u^i±Uv±, Smn = ^vm±vn±, (18) 
where a is the surface density of dust matter. For the metric (1), we have 

Ko± = b±v± + {a±/a±)v±]^±/a±, (19) 

and 

Kl± = = [<J'/'J + {a±/a±)v±]j±/a±, (20) 

where K}^_^_ = q^^KLN±, and a dot and a prime denote d/dri± and d/dx±, respectively. 

By substituting these equations (19) and (20) into the above Eqs.(13) - (15), we can obtain equations for v± 
and (7. By solving these equations, Sakai ct al.[3] derived their time evolution and showed that \v±\'^ w 10~^ <C 1. 
Moreover the following constraint equation for a is obtained: 

[7(0-' + vHR)]^ = -^K^Ra, (21) 

where R = a+a+ = a-U-. 

In the following the suffice ± is omitted except for the case when it is necessary. Here let us show the relations 
necessary for the later calculations: 

900 = -1, 922 = 933/ sin^ S = {aaf, (22) 
gOO = -l, q^^ = q^^ sinH = {aa)-\ (23) 

Koo = -{-f/a)[-f^v+ia/a)v], K22 = K33/ sin^ 9 = a-fa^[a' /a + {a/a)v], (24) 
where K = q^^ Kmn = {-i/a)[2a'/a + 3(d/a)t; + -f^v] and = Kqo, K'^^ = K^^ sin^ 9 = K^iliaaf. 
Koo - Kqoo = 2{'-f / a)[(j' / a + {d/a)v], 

K22-Kq22 = [K33-Kq33]/sm^9 = -a'ya^[a'/a + 2{d/a)v + j^v]. (25) 
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III. PERTURBATIONS AND THE JUNCTION CONDITIONS 

The perturbations of g^v and Zi^ around the background (^^jy and Z'^) are considered: 

5/.. = 9^.. + Sg,.. , Z^ = Z^^ + 8Z>'. (26) 
The perturbed tangent vectors are expressed as 

^(M) — ^(M) ~ ^^ze-^M)- (^''') 
For the perturbations of the induced metric: 

IMN — Qmn + 5qMN , (28) 

Mukohyama derived the expression 

Sqmn = e-lMflN){^9iJ.v + Cszgixv)- (29) 

The perturbed normal is expressed as 

5n^ = ^nf^n'^n^dg^x + e^^^g^^n^A^e'^^) (30) 

and 

SKmn = ^e'^Mfwi'^^szguv + CszCnQ^v - 2n^5Txij.u], (31) 

where 

= ^{Sgxn;,^ + 5g\v-^ - 5g^^.x). (32) 



This perturbed extrinsic curvatures are reduced to 



SKmn = -^IKmn - -jn^^{M)'^'(N)J>^tJ.v, (33) 



where 

I = n''n"(%. + 26 Z^,,^) (34) 

and 

JxiJ.v = 25Txij,v + SZx;ij.v + SZx-vij. + {RaiJ.\v + Rav\iJ,)SZ°' . (35) 

Under the gauge transformation 

xf" x'^ = a;'* + (36) 

5 and SZ^ transform as 

^ - - SZ^^dZ^' + e, (37) 

and it is found that Sqmn and SKmn are invariant under this gauge transformation. The junction conditions 
for the perturbations of metric and extrinsic curvature which was shown by Mukohyama[15] are expressed as 

SqMN+ = Sqmn- (38) 

and 

[SKmn]^ = -k\SSmn - qmn), (39) 
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where 



and 



5k MN = SKmn - ^ {K^SqLN + K^Sqem), (40) 

SSmn = SSmn - -^i^M^QLN + S^dqLM) (41) 
with S = Smn and 5S = q^^ SSmn- The equation (39) is consistent with a relation 

[SKmn]^ = -k'^{SSmn - ^{SSqMN + SSqMN)}, (42) 

so that 

- k^SSmn = [SKmn - SKqMN - KdqMN]^, (43) 

where SK = Sq^^KMN + q^'^'SKuN = q^^'SKuN - K^^'SquN- 

When the metric perturbations in the two regions V+ and V~ are given, we can determine SqMN, SKmn and 
SSmn using Eqs. (29), (31) and (43), and find the conditions which are imposed on the shell displacements, 
using Eq.(38). 



IV. PERTURBED QUANTITIES IN THE BOUNDARY SHELL 

In this section we express first the perturbations of the induced metric and the extrinsic curvature, using the 

gauge-invariant variables representing metric perturbations and the shell displacement, and derive the pertur- 
bations of the energy-momentum tensors in the shell. The metric perturbations and the harmonic functions in 
the homogeneous regions are shown in Appendix A. 

A. Metric perturbations 

The metric perturbations in Eq. (29) are rewritten in terms of x'^ as 

SqMNdy^Uy^ = {Sg^, + SZ^,, + 5Z,,^)dx''dx\ (44) 

where SZfj, represents the displacement of the boundary shell. SZq and 6Zi arc here treated as spatially 3- 
dimensional quantities in accord with the metric perturbations. In principle they are the quantities given in the 
shell. In Mukohyama's formalism[15] the displacements of the shell were regarded as functions of only intrinsic 
variables r, 9 and ip which are spatially 2-dimensional. In order to recover the local nature of SZq and SZi in 
our case, we impose later a localization condition, imder which their behavior is constrained so that their values 
may be nonzero only in the neighborhood of the shell. The coefficients in the right-hand side of this equation 
are expressed as follows using the gauge-invariant variables. 



1. Scalar perturbations 



Sgoo + 2SZo;o = j k^dk[-2a^^A + 2a{(t)o/ay]Q, 
Sgoi + SZi.o + SZo-i = j k'^dk[a'^{^L/a^y +MQ,i, 

Sgij + SZi.j+SZj.i= J k'^dk[{2a'^^H -2^(t)oHjQ-k(l>L{Qiij + Qjii)], (45) 

where the suffix \i denotes the covariant derivative in the 3-dimensional space with dl"^ = "jijdx^dx^ , SZo and 
SZi are expanded in terms of harmonic functions as SZo = J k'^dk zoQ and SZi = J k'^dk z^Qi, and k is the 
wave-vector, whose length is A; = |k|. Here, ^a,^h,<Po and ^i, are the gauge-invariant variables defined by 

^A = a + Ib + 1^b~^(ht + -h, 

k k a 
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(f>o = zo + iayk){B-HT/k), 
(t>L = -ZL/k + a^Hr/k^, 



(46) 



where various expressions for the metric perturbations and harmonic functions are shown in Appendix A. Here 
we use normahzed harmonic functions. 

Since dx° = Fdy'^, dx^ = Fvdy^, dx"^ = dy'^ and dx^ = dy^ with F = e°Qj = 7/a, we obtain 

^500 = J k^dk F^iCoQ + 2vClQ,i + v^{Cll1iiQ + 2fc'</>LQii)], 

SqOa = j k'^dk F[aQ,a + '^vk^(l>LQla], 

Sqab = j k^dk[QLLlabQ + '2k^HQab], (47) 

where Co = -2a^$A + 2a(^o/o)', Cl = o?{<i>L/a?)' + </>o and Cll = ^q^^h - 2{d/a)<pQ - ^k'^(j)L- For I defined 
in Eq.(34), we obtain 



2a^$H - 2-(po 



1% 
a 



{a/jfl = J k'^dk 
Equations (47) are rewritten using spherical harmonics as follows: 
5qoo = j k^dkF^ | [CoH + 2vCL'n' + v 



Q + 2(j>LQM + 2v[a^{^L/a'^y + MQa + v'^[-2a'^^A + a{4>o/ay]Q^ . 

(48) 



CllU + 2ci>L ( n" + y n 



Sqab = k'^dk 



LUYab + {ClL + 2(t)L 



a 



2a2 



n 



labYlr, 



(49) 



where Q = Iii{k,x)Yim{0,'^), the suffix / in 11; is omitted, Yim{0,^) is the spherical harmonics and Yab is a 
2-dimensional traceless tensor defined by 



^ab — ^lm\[ab ~t~ 2(J^ Tafe^m- 



(50) 



Here a and h take the value 2 or 3 and || denotes the covariant derivative with respect to the space with 

dVl^ = de^ + sin^ edip^ . 

In these equations, dots mean d/dr]+ and d/drj- in V+ and V^, respectively, primes mean d/dx+ and d/dx- 
in V+ and V~, respectively, n+ and H_ are equal to n(x+) and H(x+), and a suffix A; is omitted here and in 
the following. Here, if we change the time variables r]+ and r/- to the common variable r(= y°) in the shell by 
d/dr]± = {a±/'y±)d/dT, then we have 



Co = -2a2$^ + -a2(0o/a),r, 
7 

Cl = a^i(j)L/a^),T + 4>o, 
Cll = 2a'<^H-2a,r(l>o-^k''(t^L- 



(51) 



2. Vector perturbations 

The coefficients of Eq.(44) are 



6goQ + 25Zo-o = 0, 

5goi + SZi.o + 5Zo-i = j k'^dk[-aH + a^{(f>T/a'^y]Vi, 

Sgij + 5Zi.j+5Zj.i = J edk[-2k(l>T]Vij, (52) 



where SZq = and SZi = J k'^dhzrVi, and ^ and cpT are the gauge-invariant variables defined as follows: 

1 .2 Ml) 



k 



(53) 



The induced metric perturbations are 

Sqoo = -2F\ j k'^dk{[-aH + a^{(t)T/a'^y]Vi+vk(l)TVii}, 



Sqoa = F j edk{[^aH + a'{(t)T/a') ]Va-2vk(l>TVia}. 

Sqab = -2 / k^dk {k(t>T)Vab, 



where Via = -M'^i.a + R^{Val R^),i]. For I we obtain 



{ahfl = -^j k^dk{v[aH + {^t / a'^y]Vi + k^rVn}. 



(54) 



(55) 



Equations (54) are rewritten using spherical harmonics in the cases of odd and even parities as follows: 

(2-1) The odd-parity case 
Since Vi = (cf. Appendix A), we have 

Sqoo = 0, 

6qoa = fJ k^dk{[-aH + {a^M{cPT/a^),r]Va+vMVa,i-2{(T'/a)Va]}, 

Sqab = j k'^dk(l)T[Va\\b + Vb[\a], 



(56) 



where (1^2,^3) = (T(x)n(— Yj„i^3/sin^, y;TO^2sin^). These 3-dimensional vector perturbations are also 2- 

dimcnsional vector perturbations. 



(2-2) The even-parity case 
Harmonic functions are 



Vi = 1(1 + 1)01/ a)Yi^, 

{V2,V3) = iUay{Yim,2,Yim,3), 



(57) 



and the intrinsic metric perturbations are 



5q( 



'00 



= -2F\ j 



k^dk 



IT _ 
IT ~ 7 



Vi, 



-a^* (aV7)(^T/a^),r - V(j)T 
6qoa = fJ k^dk{[-aH + (aV7)('^T/a2),r]K + v<Pt - 2{a'/a)Va + , 

Sqab = J k^dk(f>T[Va\lb + Vblla + 2jab^Vi]. 

These expressions are reduced to 



(58) 



Sqo 



= -2F\ j 



k^dk 



n' a' 

-o^* + (aV7)(0T/a^),r -v<j>T[-^- — 



l{l + l)in/a)Yim, 



Sqoa = F j k^dk{[-a'^ + {ay^){<i>T/a'),r\{nay +vM{^'yy' -'^{'^'l'^){^'^y + 

Sqab = J fc'rf/e.^T{2(na)'i;6 + ^^^^[20^' - {Ilay]}jabYim. (59) 
In this form the dependence on spherical harmonics is found to be quite the same as that of scalar perturbations. 
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3. Tensor perturbations 

The non-zero coefficient of Eq.(44) is 

s,,.sz.,,.sz,,. = -,^je,tHipa,. (60) 

and the induced metric components are 

5qoa = -2a'^Fv j k'^dkH^^da, 

Sqab = -2a2 J k^dkH^^Gab. (61) 

For / we obtain 

{a/-ffl = -2a^ J k^dkH^^Gii. (62) 

Equations (61) are rewritten using spherical harmonics in the cases of odd and even parities as foUows: 

(3-1) The odd-parity case 
If we define a 2-dimensional axial vector Wa by 



for a = 2, 3, we obtain 



Sqao 
Sqoa 

Sqab 



Wa = {-Yim,3/ sin e, Yirn,2 sin e) (63) 
0, 

Fv J k^dkHl^\l - + 2)mVa, 

j kHkH^^\lla^)'{W,\\^ + W^^). (64) 



When we compare Wa with Va in the odd-parity case of vector perturbations, we have the relation Va = {Ji(j)Wa-, 
so that with respect to the two-dimensional angular dependence the tensor perturbations in the odd-parity case 
have the same form as vector perturbations in the odd-parity case. 

(3-2) The even-parity case 

5qoo = {Fvf j k''dkH^^L{Il/a'')Yim, 

Sqoa = k''dkH^^\l-l){l + 2)^^Yim,a, 



Sq^f, = j k^dkH^^ 



2grYab - ^njabYln 



(65) 



where the definitions of L and ^" are given in Appendix A. In this form the dependence on spherical harmonics 
is found to be quite the same as that of scalar perturbations, as well as that of vector perturbations in the 
even-parity case. 

B. Perturbations of the extrinsic curvature and the energy-momentum tensor 

in the shell 

Now let us calculate the perturbed extrinsic curvature 6Kmn using Eq.(33) and derive the perturbed energy- 
momentum tensor SSmn from Eq.(43). The former expression is rewritten as 



SK 



MN 



^7 Kmn j k'^dkJ^MN, (66) 
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where Tmn is defined by J^mn = '^^6(M)^('jv)"^^i"''- '^'^^y have the following components: 
^00 = (7/n)'[^ioo + t^(Jooo + 2Jioi)+i'^(Jiii+2Jooi)+w^Joii], 

^22 = Ji22-^vJq22-, 

•^33 = -^133+^^-^0335 

■^23 = Jl2S, + vJii23,, 

^02 = (7/a)[-^i02 + v(-/oo2 + -/112) + J012], 

^03 = (7/a)ki03 + v(Joo3 + -/113) + t^^-/oi3], (67) 

where ^20 = -^02,^30 = -^TjS) and T32 = ^23- The expressions of Ja^i/ in terms of gauge- invariant variables of 
metric perturbations and displacements of the boundary shell are shown in Appendix B. 
The energy-momentum tensor in the shell is expressed as 

— k^SSmn — [©mat]''', 

Qmn = j k'^dk^MN = ^I{Kmn - KqMN) ^ \~ j k'^dk^MN - KSqmn, (68) 

where 

^MN = ^MN — Q^^^KL Qmn (69) 

and their components are 

^00 = 7— T2 [-^122 + J022 + (-/133 + V J033)/ sin^ 6] , 
{acr) 

T22 = {<Jl?[Jwa + v{Joaa + 2Jioi) + + 2Jooi) + v/^/on] - (Ji;,3 + vJo33)/sin2 9, 

^33 = (0-7)^ sin^ 6i[Jioo + v{Jooo + 2Jioi) + v^{Jiii + 2Jooi) + v^Jqu] 

-{Ji22 + vJo22)sin^e. (70) 

The components of ©mat in terms of the gauge-invariant variables are obtained through lengthy calculations 
and are shown in Appendix C in the expanded form with respect to the shell velocity v. 



C. Junction of perturbations 

From now let us discuss the junction of metric perturbations given by [Sqmn]^ = in Eq.(38), in order 

to analyze the displacements of the shell. If the displacements are obtained for given metric perturbations in 
both regions, we can determine the perturbed energy-momentum tensor in the surface {SSmn) using the above 
equations (68). 

The following analyses in the shell are performed separately in the odd-parity and even-parity cases of metric 
perturbations, and all perturbations in each case are treated together in spite of the difference in their 3- 
dimensional behaviors. This is because the 2-dimcnsional angular dependence of metric perturbations in each 
case is the same and indistinguishable, as was shown in the previous subsection. 



Odd-parity perturbations 



Here we consider vector and tensor perturbations together, and from the condition [(^^mat) (vector) + 
((^f/MAr) (tensor)]^ = 0, wc obtaiu the following two independent equations: 



j k^dk {[-7a* + a'^{(l)T/a^),r]01cT) + ^v[(l)T{U'a - Ha') + {I - 1){1 + 2)H^^n]^ 
j k^dk^(l)T<jH + H^\a^U)'^ 



= 0, 



(71) 



where the ./;"(= r]) derivative was changed to the r derivative by d/drj = {a/j)d/dT. 

(2) 

The metric perturbations ^ and H!^ are given as functions of r and k in both regions. In order to recover 
the original local nature of ^t, on the other hand, we impose the following localization condition on it using a 
function d{x)- 



I 



(f)T±H.±a±k'^dk = (f)T±{T)d±{x), 



(72) 
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respectively, where (f>T± (t) arc functions of only r and the functions d± (x) corresponding to the shell boundary 
(X = Xb) S'lid the small width e arc defined as 



d+{x) = ^ for Xb<X<Xb + e and for < x < Xb or x > Xb + e, 



(73) 



and 



c^-(x) = l for Xb-e<X<Xb and for < x < Xb - e or x > Xb- (74) 

The above condition means that the total displacement (/ (j)T±^±cr±k^dk) has the nonzero value only in the 
neighborhood of the shell. Here we assume that e is small enough compared with X5- Then we find that (^t± 
are inversely expressed as 



- 1 r°° 

(pT±=(pT±{T)-j-^ n±(x)CT±(x)rf±(x)rfx - (pT±{T)ek-^n±{k,Xb) 
in terms of ^t±{t), using the orthonormal relations ([21, 22]) of 11 which are expressed as 

dxn±(fc,x)n±(fc,x)a2(x) = <5(fc-fc), 
dkiL±{k,x)ii±{k,xMxMx) = S{x-x) 



(75) 



Jo 



(76) 



for n± with common I. Here it should be noticed that the orthonormal relations require the integration interval 
[0, oo] of X) but it is formal and independent of the physical interval of x in the two homogeneous regions. 
Eq.(75) shows the k dependence of 0t±- Under the above condition, the two unknown variables 0t+ and 4>t- 
are determined by solving the above two equations, as follows. 
First we obtain from the latter equation of Eq.(71) 



= 0. 



(pT- =<pT++ J k''dk[H^^\a^uy]^. 
Next let us differentiate the latter equation with respect to r. Then we obtain 

k^dk [(j)T,r{Ii(T) + H^liUa^y + [^TiUa]' + H^\Ua^y]-fv/a} 
Eliminating 0T,r from the former of Eq.(71) and Eq.(78), we obtain 

J k'^dk \24>T [^(T + l{v/a)a')\ H + ^a^Ua + //^^^(Hct^)' + ^{v/a)H^^\(Jlay -{I- + 2)H]} 

Substituting Eqs.(72) and (75) to this equation, moreover, we obtain 
20T (^CT + 7(t;/a)£)] ^ + j k'^dk {la^Ha + H^^l{\la'^)' + ^{v / a)Hl^\{Il(T'^)" -{I- + 2)n]} 



where 



£ = y dfc y dxdixMxMx) n(x)^' 



(x). 



(77) 
(78) 

= 0. 
(79) 

= 0, 
(80) 

(81) 



In Appendix D, C is rewritten and reduced to a compact form. Eliminating using Eq.(77), therefore, we 
obtain 



2(j)T+ 



edkH^\lla') 



r 1 

/ k^dk {laMcT + H!^l{ILa^y + j{v/a)H^\{ILay -(/-!)(/ + 2)n]} 



(82) 



Thus Eqs.(77) and (82) give us (f>T+ and (f>T- in terms of given 3-dimensional metric perturbations. Here let 
us assume that v vanishes for simplicity. This condition holds in a good approximation, because v (w 10~^) is 
very small in our background models, as was shown by Sakai et al. (cf. §2). Then (px is reduced to 



(f)T+ = 





-:{ 


j k^dk a^Ha 


± 

+ 


. a . 









I 



k'^dkH^'lilla^) 



k'^dkH^^'illa^y 



(83) 
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2. Even-parity perturbations 



Here we consider scalar, vector and tensor perturbations together, and from the condition [(5(j'MAr)(scaiar) 
(^9mjv) (vector) + (^^mat) (tensor)]^ = 0, wc obtain the following four independent equations: 

F [ k'^dklan + 2v<f>L[^' - -n] + [-a'* + (oV7)(0T/a2),,](n<7)' + i;0T[(n<7)" - 2-{Uay 
J ^ (T a 

+1(1 + l)(n/a)] + vH^\l - + 2)(n<7)7CT} 



= 0, 



j k^dk[CLL'n + 2(j>L 



0, 



k^dk\^<l)Lll + 2(j)T{nay + 2giH!j^'>^ 

3 



n'- 



= 0, 

^(1 + 1) 

2a2 



n 



T^[(n.)'-2a'n]-^Ffn} 



= 0,(84) 



where the first and second equations come from the (00) and (Oa) components of the condition and the last 

two equations come from (ab) component. In these equations the metric perturbations <^>a, '^h, ^ and Hj, are 
given as functions of r and k at both regions. 

Now in the same way as in the odd-parity case, we impose a localization condition for (po^tl^L and 4>t as follows, 
using the function d±{x) (defined in Eqs.(73) and (74), so that the original local property of the displacement 
may be recovered: 

j <t)a±Ii±k^dk = (t>o±{T)d±{x), 
j ct>L±Il±edk = 0L±(r)d±(x), 
bT±{Il±a±yk^dk = ^T±{r)d±{x), 



I 



(85) 



where 0o±('''), ^l±(t) and 4't±{t) are functions of only t. This condition means similarly that the total 
displacements have the nonzero values only in the neighborhood of the shell. Inversely 0o± and (j)L± are 
expressed as 



k± = (t>o±{T)^ / n±a-\d±{x)dx ■^^ (t>o±{T)ek '^'n.±{k,Xb)o±{k,Xb), 



<I>L± 



h±{r)-^ U±<7ld±{x)dx^h±{T)ek-^n±{k,Xb)cTl{k,Xb)- 



(86) 



Under this condition wo find that for the above four equations there arc six unknown variables (/)o± , (t)L± and 4>t± , 
representing the displacement of the shell. But they are redundant, because only two kinds of displacements 
are allowed in the even-parity case of the 2-dimensional space, as in Mukohyama's treatment. Here SZa from 
4-dimensional scalar and vector perturbations give 4>L^Yirn.a and 4'T(JicryYim.a, respectively, in the similar 
manner. Accordingly we impose an additional condition to relate 4>l± and 4>t± respectively in the form of 



j (t)T±{^±(J±yk'^dk = j (j)L±a±k^dk or (f>T±{T) = H±{'r)- 



(87) 



Then we can eliminate ^t± using Eq.(87) and determine the remaining four variables solving the above four 
equations as follows. 

Here let us assume t; = for simplicity. Then the above equations lead to 



dkk'^^A^+{(i>o/a),-, 



0, 
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2(^L,, - 2^^l) + -4,o-a [ dkkHilLa)' 

a a J 



0, 



3^1, + 2 / k'^dkHl^^gi 
L 



= 0, 



where 



i{i+mix) 

(n(x)a(x))' 



The obtained solutions for (j)L± and 4>o± are expressed as 
= k'^dklH^^Qi]^ + H+, 



k^dk 



2 /ay\ 

3 V a /- 



and 



?!'o+ = 



yo- = 



7-2 - 4^a+,r], 



where 



A 
B 



ijdrj k'^dkl^All]^, 

I 



k^dk 



4 ^ 



(88) 



(89) 



(90) 



(91) 



k''dk[H}^>gi]^ + ^L+\M-, 



(92) 



where the definition of A4 is shown in Eq.(89) and A4± are estimated in Appendix D. Thus we can determine 
the displacements of the shell when the metric perturbations in both regions are given. 

As above-mentioned, wc can derive the energy-momentum tensor in the shell by substituting the obtained 
(t>Q,4>L and <t)T to Eq.(68), in which (j)o,(j)L and ipT can be replaced by 0O) and (px, and we find that the results 
are independent of the value of e, as long as e/xb << 1- 



V. CONCLUDING REMARKS 



We have studied the junction conditions imposed on the metric perturbations in the shell boundary in a 
four-dimensional model with two homogeneous regions . As the result we found that the perturbations with 
different parities can be separated and treated independently, and that due to the inhomogeneity in the shell 
the mode coupling is caused among perturbations with various types and different wave-numbers in each parity. 
The displacements of the shell and the perturbed energy and momentum density for given metric perturbations 
in both regions can however be consistently derived. 

Our results do not mean that arbitrary perturbations in both regions can always be adjusted, because it 
depends on the physical situation of the shell whether the obtained energy-momentum tensor in the shell is 
allowed or not. For instance, if we consider an expanding or collapsing star with the empty external region, 
any energy and momentum are not expected to be stored in the boundary and only gravitational radiation is 
released in the external region, because there are no matter perturbations. 

In the present study, though we treated the displacement of the boundary shell as spatially 3-dimensional 
quantities, wc could recover the original local property of the shell displacement (that it has the values only 
in the neighborhood of the shell) by imposing the localization condition and could solve the equations for the 
displacements. 

Our result will be useful to the analyses of CMB anisotropy observed in the nonlinear structures such as 
overdense regions (such as superclusters) and underdense regions (such as voids) in which we may live. 
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APPENDIX A: PERTURBATIONS IN A HOMOGENEOUS MODEL 

We show here the expressions for metric perturbations and harmonic functions in a homogeneous model with 
the line-element in the text, paying attention to the 2-dimensional parity operation P : 9 ^ it — 6, (fi ^ n + ip, 
where x"^ = 9 and = <^.[13, 20] 

1. Scalar perturbations 

Metric perturbations are expressed as 

500 = -a^[l + 2A{r])Q], 
goi = -a'^B{r])Qi, 

gij = a^{[l + 2HL{v)QHJ+2HT{v)Q^J}■ (Al) 

Here Q is a scalar harmonic function with the wave-number k satisfying the equation 

(A + k'')Q = 0, (A2) 

where the Laplacian A is defined by Acp = (j)\i(t>\j and \i denotes the three-dimensional covariant derivative 
in the space dJ? = 'fijdx^dx^. Qi and Qij are defined by 

Qi — k Q\i, 

Qij = k-^Qiij + -jijQ. (A3) 
Harmonic function Q is expanded using the usual spherical harmonics Yim as 

g = nr(x)i;m(0,^), (A4) 

where k"^ = 'n? — K for the spatial curvature K = 1,0, —1. Then from Eq.(A2) we obtain 

H" + 2{a'/a)n' + [n^ _ l{l + l)/a'^]Il = 0, (A5) 

where we omitted siffices of H" for simplicity, a prime denotes the partial derivative with respect to a;^ = % and 
cr(x) = sinx,X)Sinhx for K = 1,0,-1, respectively. The normalized expressions of 11 are found in Wilson's 
paper[21]. For the above-mentioned parity operation, we have 

P{Yi^) = {-l)%m. (A6) 

This property is called the even parity. 

2. Vector perturbations 

Metric perturbations are expressed as 

500 = -a^, 

501 = -a2B(i)(r;)y„ 

Qij = a''[^ij + 2H!^\'n)Vi,]. (A7) 
Here is a vector harmonic function with the wave-number k satisfying the equations 

(A + fc2)v;=0, V; = 0, (A8) 

and Vij is defined by 

Vij = -{2k)-\V,\,+V,\i]. (A9) 
In vector perturbations there are two cases with different parities. 
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a. The case of odd parity 
The components of vector harmonics Vj are expressed as 



Vi = 0, 

V2 = -a nYim,3/ sm9, 

V3 = anYim,2sme. (AlO) 
In this case we have P{Vi) = (— l)'~^Vj and the property is caUed the odd parity. 

b. The case of even parity 
The vector harmonics Vi in another case is 

(1/2,1/3) = (na),i(>';™,2,i1™,3)- (All) 

In this harmonics, we have P{Vi) = (—1)'!^, so that the property is called the even parity. 

3. Tensor perturbations 

Tensor harmonics Gij satisfies 

G.fcj^: = -fc'G.fe, „ = 0, Gl = 0, (A12) 

and the components are expressed as follows using 11" and Yim- The tensor perturbations also have two cases 
with different parities. 

a. The case of odd parity 



Gil 


= 0, 




G22 


= {Ua 


^)A-Xim/sinO), 


G33 


= {Ua 


^),iXim sin 6*, 


G23 


= (Her' 


^),iZijn sin 6*, 


G12 


= (l- 


i)(; + 2)n(-yj„,3/sin 


Gi3 


= (l- 


l)(Z + 2)nii„,2sin^, 



(A13) 
where 

Xirn = 2(r,„,23 - cot 0Yim,3), (A14) 
Zlm = Yim,22 - COt 6Yim,2 - ^^m.Ss/ sin^ ^- (A15) 

If we define a two-dimensional vector Wa with the components W2 — —Yimfij smQ and W3 = yim,2sin0, it is 
found that Gab and Gio is reduced to 

G„6 = (<72n),i(w-„||(, + W6||„), 

Gi„ = (Z-l)(Z + 2)nW„. (AI6) 
For the parity operation, we have the odd-parity property P{Gij) = (— l)'~^Gij. 
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b. The case of even parity 

The tensor harmonics in another case are 

Gil = — o-nYim, 
(G22,G33) = -^IlYirn{l,sm^e)+grZirn{l,-sm^0), 

G23 = Ql^im, 

G12 = (/-i)a + 2)i(ncT),ii,™,2, 

a 

Gi3 = {l-l){l + 2)-{\la),^Yirn,s, (A17) 
where L = l{l + - + 2), and 

= aa'n 1 + [i + Z + 2) - _ 2)a2]n. (A18) 
In this case we have the even-parity property P{Gij) = (— l)'Gij. 

APPENDIX B: Jx^,^ 

The resulting expressions for J\^i, for three types of metric perturbations are shown here. 

1. Scalar perturbations 



Jooo = 2{-a^^A + 4>o-S^o + m^f -^]MQ, 

Jioo = 2[aHA + h-^lk + ['2{^f ~^]MQ,i, (« = 1,2,3) 

Jioi = 2{a2$H-[-0o + (--3(-)2)</,o]}Q + 2(<^i-4-</,i)g,n, 

Ct Ct Qj CI 

JiOb = 2{^L-2-^L)[Q.ib--Q.b], (& = 2,3) 
a ^ . . 

Jill = [2a2$H-4-0o-2-(</)i-2-<^i)]Q,i+2(/.LQ,iii, 
a a a ^ , 

J116 = (2o2$H-2-(/)o)g,6 + 20i[g,n-2-Q,i+2(-)2Q],6, (6 = 2,3) 
a (J (J 

J122 = -2{^Ha^ + -hyQ,i + 2(Ai{Q,i22 - 2-Q,22 + cra'g.n - [{era')' - 2(-)V2]Q,i}, 
a cr a 

a g' 
J133 = -2(4>ffa2 + sin^ 610,1 + 2</)l{Q,i33 - 2— g,33 + sin^ WQ.n 

a (T 

— sin 0[(^a')'-2(-)2a2]g,i + sin 9 cos 0(Q,i Q),2}, 

J123 = -2</.Lcr2 sin 6'[g/(a2 sin 0)],i23, 

Jon = [-2(a2$^)-+4ad$A-2-(/)o]g + 20o{Q.ii + [- + (-)']Q}-4-(/)LO.ii. 

a a a ^ a 

J022 = [-2(a2$jj)+4aa$A-2^<^ok2Q^2<^o{Q,22 + [^ + (^)']a2g + aa'0,i} 

- 4-0z,[g,22+(7a'O,i], 

Jo33 = [-2(a2$H)+4ad$A-2-0ok^sin2^g + 2<^o{g,33 +[- + (-) V'sin^^Q 

a . a a 

+ ctct' sin^ ^Q,i + sin 61 cos 61(3,2} - 4^(/>z,[Q,33 + era' sin^ 6ig,i + sin0cos^g,2], 

J023 = 2(00 -2^(^L)(g,23- cot 0Q,3), 



h a' 

Joib = 2(0o-2-(/.L)(Q,ib Q,b), (6 = 2,3) 

a a _ 

Jooi = [-2aHA + 2{^o-2^(l>o)-2^{^L-2^(l>L)]Q,i, (i = 1,2,3). 

2. Vector perturbations 



Jooo = 0, 

Jwo = 2[4>T-6-4>T + 2{-- + 2{-f)(j>T-a'{i> + -^)]Vi, 
q a a a 

Jioi = 2{4)T - 2^(I)t)Vi,i, 

Jwb = a2*(H,i - Vi,b) + 2(<^T - 2-<f>T){Vi,b - -Vb), {b = 2,3) 

a a 

Jill = [2ad^ -2-{0t-2-<Pt)]Vi+2^tVi,ii, 

a a ^ 

Ji22 = 2-a^iaH - (Pt)Vi + 2<^t[Vi,22 + <7a'Vi,i - 2— ¥2,2 - {cTaJVi + 2(-) Wi], 
a a a 

a ■ cr' 

J133 = 2-a^ sin^ e{aH - ^t)Vi + 2^t[Vi,33 + (tct' sin^ eVi,i - 2— ¥3,3 
a ' ' ^ a ' 

- (aa'y sin^ 61^1 + 2{- fa^ sin^ eVi + sin6'cos6'(yi 2 - 2—V2)], 
^ a ^ ' (7 

J123 = 2(t)T[Vl,23-—{V2,3 + V3,2) + COte{-Vi,3 + 2 — V3)], 

Jiib = 2(l)T[Vi,ia--{Vb,i+Vi.b) + 2{-fVb], (& = 2,3) 
cr ^ (T 

Jooi = [2aa^ -2'^{^T-2%l>T)]Vi, (i = 1,2,3) 
Jon = -2{aH + 2^(l}T)Vi,u 

J022 = -2(o^* + 2^0t)(V2,2 + (T(t'Fi), 

J033 = -2iaH + 2%T){V3,3 + crcr'sm^eVi), 

J023 = -{a^'i + 2-cl)T)iV2.3 + V3,2~2c0teV3), 

a ^ 

J016 = -{a''^ + 2-(j)T){Vi,b + Vb,i-2-Vb). (6 = 2,3) 
a a 



3. Tensor perturbations 



Jooo = Jqoi = J002 = Joo3 = JlOO = 0, 

Jio,- = -2{c?iif^)Gij, (i,i = 1,2,3) 

Jill = 2a^ii^^^ G\i,\^ 

J122 = 2c?ii^^\2Gi2^2-G2i,x+2aa'Gii\, 

J133 = 2a2ij^^^[2G'i3,3-G33,i + 2cra'sin2 6'Gii+2sin6'cos6'Gi2], 

J123 = 2a2jr^'^[Gi2,3 + G3i,2-G23,i-2coteGi3], 

J116 = 2a2ij^'^[Gii,b-2-Gib]. (6 = 2,3) 
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APPENDIX C: THE INTRINSIC QUANTITY Qmn IN THE SHELL 

The expressions of ©mat are shown for three types of metric perturbations, where M and A'' take the values 
0, 2 and 3. 

1. Scalar perturbations 



where 



e 



00 



a 



2 ( <t' 1 2 (t' 

-a\<pH + 2$a)Q + a'$HQ,i > + ^-(-200 



2 d • 



1 



1 



a 



+ --0lQ,i - ^0l{^[(Q,i - 2-Q),22 + ((3,1 - 2-Q),33/sin^ Q\ 
- 2[— + - 2(-)2]Q,i + cot^(Q,i - 2-Q),2} 



^0 



a' 
2 d • 

+ ^-(i-o 



a 



/(/ + 1) 



a a 
2v , a 



2-+4(-)2-^^)Q + 4^Q,i 



2 d , • a , 

A^-v + v){(po (Pq)Q 

a"^ a a 



-^H-[Qm + ^(Q,22 + Q,33/sin^^) 



9 



+ 4— Q,i + ^cot(9Q,2] - ^— +0(t;2), 
a a-^ a'* cr 

1 .ct' 



22 



e, + -(/.l{(Q,i33 - 2-Q,33)/ sin^ 6 + cot 0(Q.i2 - 2-Q.2) - 4-Q,22} + 0(f), 
a a a a 



633 = e, sin2 e + -c^l{{Q.122 - 2-Q,22) sin2 e - 4-(Q,33 + sin0cos0g,2)} + Q{v), 

Oob = 2i^L-2-(l)L){Q,l + —Q),b-^ — (l>oQ,b-^v{^L-'2-(l)L + MQ,b 

a a (7 a a a a 

^1 



V 



. a 



a 



, a 



-^{[2a\^H - ^a) + 2(/)o + ^-{H - 2-(/)i)]Q,6 + 2</)i[Q,n + 2-Q,i 



a 



a 



a 



- 2{-fQlb} + Q{v% (6 = 2,3) 

cr ^ / • 

©23 = -WL[Q.i2-2-Q.2-cot0(Q,i-2-Q)].3-^;(0o-2-</.L)(Q,2-cot0g),3 
a ^ cr cr a 

+ 2(2- - 3-t; - v)<^l(Q,23 - cot ^©,3)} + 0(^2), 
cr a 



(CI) 



e,/a2 = -ia2(cI>^ + $^)Q.i-5-i(a2$^--(/)o)Q-;i(^i-2-<^i)Q,i 
a a a a la a. 

+ -</'L{--g,ii + [- - 3(^)2]g,i + - (^)2]Q,i} + 0(z;). 
a a a a a a 



(C2) 



2. Vector perturbations 



©00 = -4r(-(a2*-0T)T4+</'Th2-T/i,i+4(-)Vi + ^(^1,22-2-1/2,2 
a-^ la a a cr^ V cr 

+ (^1,33 -2-1/3,3)/ sin2 + cot f?(T/i,2 -2-1/2)) +2-1/1,1 - ^^^i^Vi]} 

+ ^\2-<j)TV, 1 - 2[a.^^ + a^ct>T/a^y]-V, 
a-* L o cr 

+ (a^* + 2-,^t)[^(1^2,2 + 1^3,3/ sin^ 6) + 2-Vi] 
a a cr 

+ 4^[-aH + a^,l>T/a^y]vA+0iv^), 
cr J 



1 ct' a' g' 

622 = + -<^t[(14,33- 2-^3,3)/ sin2^-4-y2,2 + cot ^(^1,2- 2-1^2)] 

a_ (T _ fj (T 

^t(3-u + v)y2,2 - -{c?^ + 2-^t)V3,3/ sin^ 61 + 0(t;^), 

Qj Ct Qj (X 

633 = sin2 ^? + -(j)T[iVi.22 - 2— ¥2,2) sin2 0-4- (^3.3 + sin 9 cos ^^2)] - -<^t(3-u + i})(1^3,3 
a ' (T a a a 

+ sin6»cos6»T/2) - -(a^* + 2-</)t)V2 2 sin^ 6* + 0(t;2) 
a a 
1 cr' cr' 

623 = (?iT[Vl,23 + —(T^2,3 + ^3,2) -cot 61(^1,3 + 2 — ^^3)] 

a a ^ a 

+ ^[a^^ + 2(- - — )H(V2,3 + 14,2 - 2coteF3), 
Za a a 

©06 = -\{la^'^{V,^,-Vi^b) + {4>T-2-<PT)iVi,b--Vb) 

2 a a 

+ + + v)[-a2* + a2(<^y/a2)-]y4 

(J o ^ 

- - -(<^T - 2-<^t)]H + MVhib - -{Vb,i + Fi,6) + 2{-fVb]} 
a a a a a 

+ [-i^;<^y(2-+3% + *)](yi,6 + H,i-2-H) + 0(t;2). (6 = 2,3) (( 
a a a a 



where 



e„ = -^^</,r[i4,i + 4-yi]-— [^T-6-^T + 2(-- + 2(-)2)0r 

a a a a a a a 

2 • / 

- + -^)]Vi + —{-{a^^ - ^t)Vi + M-Vi,i 

a a a a 

12' I 

+ 2{-fV^ - (-)'a-Vi]} - —{\2v(^T - -0t) + #t]Vi 1 + -v\c?^ 
a a a a a 

2 • f 

+ a^{(l)T/a^y]Vi} - 2— <^t(3-i; + v) — Vi - -(a^* + 2-(^tWV^ + Q{v^). (( 
a a a a a 



3. Tensor perturbations 



eoo = -a-'ir^'^[a-2(2Gi2,2-G22,i + (2Gi3,3-G33,i)/sin2 0)+6-Gii 

+ 2cot6lcr-2Gi2] +0(t;), 
622 = aH^^^ [3c7£7'Gii + (2Gi3,3 - ^33,1)/ sin^ ^ + 2 cot ^Gia] +Q{v), 



633 = aHl^^ [3aa'Gii + (2Gi2,2 - ^22,1)] sin' 6 + 0(u), 



. (7 



623 = a^?^'^[4^G23-(Gi2,3 + Gi3,2-G23,i-2cot0Gi3)]+O(t;), 
606 = -H^T^Gu + Q{v). (6 = 2,3) (( 

APPENDIX D: INTEGRALS C AND M 

The integrals £ and M. in Eqs. (81) and (89) are rewritten here using the the orthonormal relations (76) 
First we have 

= W'{x)lc7{x)] I dxdix) j dfcn(x)a(x)n(x)<7(x) 

= K(x)/^(x)]rf(x). (I 

Since d{x) = 1 in the neighborhood of the shell, we obtain 

£± = {cj'/a)±. (I 
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Next we notice that 



nn' 



n 



(no-)' a 



Wy 



(D3) 



and the factor |ncr'/(ncr)'| is <C 1, if assume that the perturbations shorter enough than the sheU size contribute 
dominantly to the integral. Then M. reduces to 



n(x)a(x) 



a'(x)\' n(x)a'(x) fa'ix) 



- 1(1 + 1) I dxd{x)a{x) I dfen(x)a(x)"^^^''^^^ 



<t3(x) 



1 ^ 



^{xVix) 



<x) 



(n(x)^(x))' 



(a')^+/(l + l)(l-0(^) 



Ha' 



^^(x), 



(D4) 



where (Ilcr'/ (Ila-)') represents the root-mean-square of Her'/ (IIo-)' in the fc-space. Accordingly, we obtain 



Ma 



— [{a')l+l{l + l)]. 



(D5) 
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